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We evaluate numerically the three-momentum-space gluon propagator in the lattice Landau 
gauge, for three-dimensional pure-SU(2) lattice gauge theory with periodic boundary conditions. 
^^' ■ Simulations are done for nine different values of the coupling /3, from P — (strong coupling) to 

/3 — 6.0 (in the scaling region), and for lattice sizes up to V = 64^. In the limit of large lattice 
volume we observe, in all cases, a gluon propagator decreasing for momenta smaller than a constant 
value pdec- From our data we estimate pdec ~ 350 MeV. The result of a gluon propagator decreasing 
in the infrared limit has a straightforward interpretation as resulting from the proximity of the 
so-called first Gribov horizon in the infrared directions. 
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I. INTRODUCTION 
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Cn , The infrared behavior of the gluon propagator in lattice Landau gauge has been the subject of several numerical 

^^ ' studies p]-jlO||. In fact, although this propagator is a non-gauge-invariant quantity, the study of its infrared behavior 
provides a powerful tool for increasing our understanding of QCD, and for gaining insight into the physics of confine- 
ment in non-abelian gauge theories (see for example pT|). In particular, the infrared behavior of the gluon propagator 

"S I can be directly related [|2[ to the behavior of the Wilson loop at large separations, and to the existence of an area 

'~~i . law. 
O On the lattice, the Landau gauge condition is imposed [lUlq] by finding a gauge transformation which brings the 

O functional Suig], defined in eq. (|j) below, to a minimum. A lattice configuration satisfying this minimizing condition 

i-G belongs to the region 51 of transverse configurations, for which the Faddeev-Popov operator is nonnegative Q-|6|. 

This region is delimited by the so-called first Gribov horizon, defined as the set of configurations for which the 
smallest non-trivial eigenvalue of the Faddeev-Popov operator is zero. (The Faddeev-Popov operator has a trivial null 
eigenvalue, corresponding to a constant eigenvector.) 

The restriction of the configuration space to the region fi implies a rigorous inequality ||l4-Eq| for the Fourier 
components of the gluon field A^ {x) . From this inequality, which is a consequence of the positiveness of the Faddeev- 
Popov operator only, it follows that the region fi is bounded by a certain ellipsoid 6. This bound implies proximity 
of the first Gribov horizon in infrared directions, and consequent suppression of the low-momentum components of 
the gauge field, a result already noted by Gribov in Reference ||l^. This bound also causes a strong suppression of 
the gluon propagator in the infrared limit (i.e. for momentum p ^ 0). In fact, Zwanziger proved |jlq,tt8[ that, in four 
dimensions and in the infinite-volume limit, the gluon propagator is less singular than p~^ in the infrared limit and 
that, very likely, it does vanish in this limit. More precisely, in Reference p8[ it was proven that, in the infinite- volume 
limit, the "gluon propagator" D{H, p) goes to zero as p ^- for almost every H . Here H should be interpreted as the 
strength of a spatially- modulated magnetic field coupled to the gluon field A^(a;), and the standard gluon propagator 
is obtained for H — 0. A similar result holds in three dimensions: one obtains [Q|8| that, in the infinite-volume 
limit, the gluon propagator must be less singular than p^^ as p — + and that, very likely, it vanishes in the infrared 
limit. 
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A gluon propagator vanishing in the infrared hmit was also found — in four dimensions — by Gribov |^^. More 
precisely, he obtained the expression p^ /(p^ + 7), where the mass scale 7^'^*^ arises when the configuration space is 
restricted to the region fi. This propagator agrees with the zeroth-order perturbative prediction p~^ at large momenta, 
but gives a null propagator at p = 0. The mass scale 7-'^''* marks the transition point between the perturbative and 
the nonperturbative regimes. A propagator that is a generalization of the one obtained by Gribov has also been 
introduced in Reference |l9|] as an Ansatz for a non-perturbative solution of the gluon Dyson-Schwinger equation 
(DSE). 

Let us notice that a gluon propagator vanishing in the infrared limit is in complete contradiction with the p~* 
singularity obtained when the gluon DSE is approximately solved in the infrared limit |O,E0[. However, a recent 
study |2l| has shown that this singularity is obtained only if the ghost contributions to the gluon DSE are neglected. 
In fact, when these contributions are included, the gluon propagator vanishes in the infrared limit [ pi| , in qualitative 
agreement with References [[l6|-[l8|. 

In this paper we present the first numerical study of the infrared behavior of the gluon propagator in three di- 
mensions.!^ As explained in References [p3|-P5| , nonabelian gauge theories in three dimensions are similar to their 
four-dimensional counterparts, and results obtained in the three-dimensional case can teach us something about the 
more realistic four-dimensional theories. Of course, the advantage of using a three-dimensional lattice is the possibility 
of simulating lattice sizes larger than those used in the four-dimensional case. This is particularly important in the 
study of the gluon propagator since Zwanziger's prediction of an infrared-suppressed gluon propagator is valid only 
in the infinite- volume limit. 

We recall that in some recent numerical studies (in four dimensions) IqJq] a sensible change in the infrared behavior 
of the gluon propagator has been observed for momenta smaller than a turn-over value pto, in agreement with the 
prediction of a gluon propagator less singular than p^^ in the infrared limit. Also, the numerical data obtained in 
Reference H have been successfully fitted by a Gribov- like formula M. Finally, in References |0-0| we have observed 

— for the SU(2) group, in four dimensions, in the strong-coupling regime and in the limit of large lattice volume 

— a gluon propagator decreasing as the momentum goes to zero. A similar result has also been obtained recently 
]lO[ for the SU{3) group in the strong-coupling regime, both in the quenched and in the unquenched case. Let us 
notice that Zwanziger's predictions |16|Jl8|| for the gluon propagator are (3 -independent: in fact, they are derived only 
from the positiveness of the Faddeev-Popov operator when the lattice Landau gauge is imposed. Thus, results in 
the strong-coupling regime (i.e. /3 ^ 0) are a valid test of these predictions. Nevertheless, it is important to extend 
this result to higher values of /3, possibly up to the scaling region. Of course, as /3 increases, one needs to consider 
larger lattice sizes in order to probe the infrared behavior of the gluon propagator. Moreover, as said above, the 
prediction we want to test applies only in the infinite- volume limit. In four dimensions it was found M-M that the 
lattice size at which an infrared-decreasing gluon propagator starts to be observed increases with the coupling. This 
makes practically unfeasible, with present computational resources, to study numerically the infrared behavior of the 
gluon propagator in the four-dimensional case and at values of (3 in the scaling region. Our hope is that in this work, 
by studying the three-dimensional case, we can consider lattice volumes that are large enough to allow a decreasing 
gluon propagator to be observed not only in the strong-coupling regime but also in the scaling region. 

II. DEFINITIONS AND NOTATIONS 

We consider a standard Wilson action for SU{2) lattice gauge theory in 3 dimensions 

3 



^ U.I/=1 X ^ ^ 



(1) 



where U^{x) G SU{2) are link variables, g is the bare coupling constant, a is the lattice spacing, and e^ is a unit 
vector in the positive /x direction. We assume periodic boundary conditions. For the matrices U E SU{2) we use the 
parametrization U = UqI- + iU ■ a, where IL is the 2x2 identity matrix, the components of it = (cr^, ct^, tr^) are the 
Pauli matrices, Uq £ ^, U £ ?ii'^ and Uq +U ■ U — 1. Notice that in eq. (nl) the lattice spacing a is necessary in order 
to make the action S[U] dimensionless; in fact, in the three-dimensional case |2^-g^, the coupling g^ has dimension 
of mass, and in order to obtain a dimensionless lattice coupling we have to set f3 = A/{ag^). 



^Preliminary results have been reported in [p2[ 



We define the gauge field A^(a;), which belongs to the SU{2) Lie algebra, as 



We also define 



M^) ^ 2 ^^^^""^ ~ ^^^"^^1 " ^U,ix)-a. 



4(a;)EE '^[A,ix)a'] ^Ul{2 



(2) 



(3) 



where a^ is a Pauli matrix. Note that a '^A^Jx) approaches {1/2) g[A^'^°'^^'^]^{x) in the continuum limit, where 
^^(cont)j6 ^^-j jg -j-j^g (unrenormalized) conventional vector potential. 

In order to fix the lattice Landau gauge we look for a local minimumg of the functional ||l|,[l3|] 

£u[9] ^ 1 - ^ E E Y [g{x)U,{x)g\x + e,)] , (4) 

/^— 1 X 

where g{x) S SU{2) are site variables, and V = N^ Nt is the lattice volume. (Here Ng is the number of lattice sites 
in the two spatial directions, and Nt is the number of lattice sites in the time direction.) If the configuration {[/^ (a;)} 
is a stationary point of the functional £u [g] then |1^] the lattice divergence of A'^ (x) is null, namely 



iV-Atix)^^. [Al{x)-Alix-e,)] ^0 



y X, b 



(5) 
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This is the lattice formulation of the usual (continuum) Landau gauge-fixing condition. By summing eq. (|5|) over the 
components x^ of x with [i ^ v^ for fixed v^ and using the periodicity of the lattice, it is easy to check R] that if the 
Landau gauge-fixing condition is satisfied then the quantities 



t{x,) ^Y.Y. ^'(^) 



(6) 



are constant, i.e. independent of Xy. 



III. GLUON PROPAGATOR ON THE LATTICE 



The lattice space-time gluon propagator is given by 

Dl^,{x~y)^{^Al{x)Al{y)). 
To go to momentum space we can use Formula (3.1a) in Reference |l6| and obtain 
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(7) 

(8) 
(9) 



Here \x goes from 1 to 3, and fc has components (fc^:, ky,kt). In our simulations we consider the values kx Ng = ky Ng = 
0, 1 and kf Nt = 0,1, . . . , Nt — 1 ; and the momentum-space gluon propagator is studied as a function of the magnitude 
of the lattice momentum 



■^Here we do not consider the problem of searching for the absolute minimum of the functional £u[g] , which defines the 
so-called minimal Landau gauge flq]. In fact, as stressed in the Introduction, the prediction of an infrared-suppressed gluon 
propagator is valid for any configuration in the region i}, i.e. for local as well absolute minima of the functional £u[g] ■ 
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If we define the momentuni-space gluon field as 

Alik) EE J2 K(^) exp [2m (k ■ x + fc^/2)] 

X 

then eqs. (g) and (g|) can be rewritten as 



Mlfc 



(10) 

(11) 

(12) 
(13) 



Notice that -D(O) in eq. (j^) [or in eq. (|l^)] is not given by !?(&) in eq. (g) [or in eq. ([l3|)] at A: = 0. The difference is 
due to the Landau gauge condition — the continuum-like condition as in eq. (m — which in momentum space reads 



Y,p,{k)Alik) = 



V k, h. 



(14) 



If k 7^ (0,0,0) we obtain that only two of the three Lorentz components of A^{k) — and therefore of A'^{x) — are 
independent. This explains the factor 6 (instead of 9) in eqs. (ph and (O). 

Let us also note that the zero three-momentum gluon propagator I?(0) can be written as 



D{Q) 



^E 



{A'f 



(15) 



where y^*^ = V ^ J^x ^m ('^) ^^ ^^^ zero-momentum component of the gluon field A'' (x) . Notice that a nonzero value 
for the constants A^^ is a lattice artifact related to the use of periodic boundary conditions and to the finiteness of 
the volume. In fact, after the Landau gauge condition is imposed, these constants are identically null — even on a 
finite lattice — if free boundary conditions are considered [Eq] , while in the periodic case they must go to zero in the 
infinite- volume limit ||l5| , [l6| . Of course, due to the volume factor in eq. (15), the latter result does not imply that 
D(0) should be zero in the infinite-volume limit. Nevertheless, as mentioned in the Introduction, it has been proven 
by Zwanziger |l(]| , |l§|| that, in this limit and in three dimensions, the gluon propagator is less singular at momentum 
p = than p^^ and that, very likely, it vanishes in the infrared limit. 



IV. NUMERICAL SIMULATIONS 

In Table | we report, for each pair (/3, V), the parameters used for the simulations .|j Overall, we have considered 
about 4100 configurations. In all our runs we have started from a randomly chosen lattice gauge configuration. To 
thcrmalize the gauge configuration {f/^(a:)} we use a hybrid overrelaxed (HOR) algorithm | 



.e. m microcanonical 

(or energy-conserving) update sweeps are done, followed by one standard local ergodic update (heat-bath sweep) of 
the lattice. In order to optimize the efficiency of the heat-bath code, we implement two different SU{2) generators 
(methods 1 and 2 described in Appendix A of Reference |Gq| , with hcutoff = 2) . In our case we did not try to find 
the best tuning for the value of ra. By analogy with the four-dimensional case [0,g| we set m = Ns/2. 
For all the pairs (/3, V), we evaluated the integrated autocorrelation timen Tint for the Wilson loops 



^Computations were performed on a SUN Ultra-1 and on a SUN Ultra-2 at the Universita di Roma "Tor Vergata", where 
part of this work has been done, and on an ALPHAstation 255 at the ZiF-Center in Bielefeld. 

''For a definition see for example [M. To evaluate the integrated autocorrelation time we use an automatic windowing 
procedure |E9| with two different window factors (6 and 15). We also employ a method pol] based on a comparison between the 
naive statistical error with a jack-knife binning error |13]. In all cases we checked that these three estimates are in agreement. 



^M ^ gV Y ^ ^ C^m;V(^) Z = 1, 2, 4, . . . , Ns/2 , (16) 

where 

f/;,;Ua^) = U^,{x) . . . U^^ix + {l-l) e^) U,{x + {l-l)e^) ... U,{x +(/-!) e,, + {I - I) e,) 

X U;\x + a - 1) e^ + (/ - 1) e,) ... U-\x +(/-!) e.) C/-i(x + (/ - 1) e,) . . . U~\x) , (17) 

and for the Polyakov loops 

AT T ^" 

^/^^-^EE Y n ^^^.(^ + ^^^6^) A^ = 1,2,3. (18) 

In all cases we obtained Tint ^ 1- (Note that Tint = 0.5 indicates that two successive configurations generated in 
the Monte Carlo simulation are independent.) Since, for all pairs (/3, V) and for all quantities, the number of sweeps 
between two consecutive configurations used for evaluating the gluon propagator (see Table |) is much larger than the 
corresponding integrated autocorrelation time, we may conclude that these configurations are essentially statistically 
independent. 



For the numerical gauge fixing we use the so-called stochastic overrclaxation algorithm | |3l| , |32[ . In all our simulations 
we stop the gauge fixing when the condition 



^Y. [iy-Af{x)\^ < 10-12 (19) 

a:, b 

is satisfied. [See eq. (P) for the definition of the lattice divergence (V • A) (x) of the gluon field ylj^(x).] This is 
equivalent |32| ] to fixing the minimizing functional £u[g] up to about one part in 10^^. In the final gauge-fixed 
configuration we also evaluate |^^ 

Q^ ^E]^E [QtM-Qt]' [qIY' . (20) 



Xtjjb 



where Q\ = N^'^ Yl,x Qti^'^)^ ^^^ Qti^i') ^^ been defined in eq. (||). The quantity Q should be zero when the 
configuration is gauge-fixed, and it is a good estimator of the quality of the gauge fixing. As in References |32[ , we 
found that the stochastic overrclaxation algorithm is very efficient in fighting critical slowing- down |29[| , and in making 
the quantities Q'l{x,y) constant, i.e. Q sa 0. In particular, by averaging over all the gauge-fixed configurations, we find 
Q — 3.5(10) X 10"^. We also obtain Q < 10~* for 69% of the gauge-fixed configurations. 

A. String tension and lattice spacing 

For each coupling f3 we evaluate the average plaquette (W^i.i) (see Table ||). Results for /3 — 5.0 and 6.0 are in 
agreement with the data reported in Table 15 of Reference [2^. In Figure Q we also plot (W^i,i) as a function of the 
coupling /3, and we compare the numerical data with the leading strong-coupling expansion /3/4 and weak-coupling 
expansion exp(— 1//3). It is clear that the crossover region from strong coupling to weak coupling occurs around 
/3 « 3, in agreement with Reference 123], and that our simulations range from the strong-coupling region up to the 
weak-coupling one. 

Following Reference [g5| we also evaluate, for /3 > 3.4, the tadpole-improved coupling j3i = (3 (VKi.i) (see Table ||). 
Then, by using the fit given in eq. (67) of that reference (which is valid for (3 > 3.0), we calculate the string tension 
■^/ct in lattice units, and the inverse lattice spacing a^^ using the input value \/a = 0.44 GeV (see the last two columns 
in Table U). 



Finally, in Table III we report (for each (3 > 3.4) the lattice spacing in fm, the largest lattice volume Vmax considered, 
the corresponding physical volume in fm^ , and the smallest non-zero momentum (in GeV) that can be considered for 
that lattice. Thus, in this work, we can explore the infrared behavior of the gluon propagator for momenta as small 
as j3 « 110 MeV, in relatively large physical volumes, and for couplings f3 above the strong-coupling region. Let us 
notice that, if we compare the data for the string tension (in lattice units) with data obtained for the SU{2) group in 
four dimensions (see for example Table 3 in Reference [p3[), then our largest value of /?, namely 6.0, corresponds to 
f3 « 2.4 in the four dimensional case. 



B. Gribov copies 

In this work we do not consider the problem of Gribov copies (see for example MM and references therein.) This 
is motivated by our finding in the study of the four-dimensional case. In fact, in References |7|g] we checked that, for 
the SU(2) group in the four-dimensional case, the influence of Gribov copies on the gluon propagator is of the order 
of magnitude of the numerical accuracy. (A similar result has also been obtained for the Coulomb gauge |3J].) In 
fact, from Table 2 in Reference ||8|, it is clear that data corresponding to the minimal Landau gauge (absolute minima 
of the functional ft/ [5]) £^re in complete agreement, within statistical errors, with those obtained in a generic Landau 
gauge (local minima of £u[g])- In particular, this seems to be the case even at small values of the couphng (3, namely 
in the strong-coupling regime, where the number of Gribov copies is higher and their effects, if present, should be 
larger and more easily detectable. 

C. Z{2) symmetry 

In References p5| it was shown for the four-dimensional case that, at very large f3, the data for the gluon propagator 
are strongly affected by the broken Z{2) symmetp^ In particular, one can consider all the possible combinations of 
signs of the average Polyakov loops (P^) [see eq. ([l8|)], for a total of 2"* = 16 different states. Then, if the expectation 
values are evaluated only over configurations belonging to the same state, the gluon propagator takes different values 
in different states |3g]. Here we did the same analysis at P — 5.0 and lattice volume V = 16'^ with 1000 configurations. 
Since we work in three dimensions, there are 2^ = 8 possible states, i.e. combinations of the signs of the average 
Polyakov loops. We obtain that, also in this case, the gluon propagator depends strongly on the state used for 
evaluating the expectation value. For example, for k = (0, 0, 1) the smallest value — D{k) — 1.50(3) — is obtained for 
the state characterized by positive Polyakov loops in the three directions, while the largest value — D{k) — 1.68(3) — 
corresponds to the state characterized by negative Polyakov loops in the three directions. The two values clearly differ 
by several standard deviations. A similar result is obtained when other momenta k are considered. This observation 
may explain why the data for the gluon propagator are usually characterized by large statistical ffuctuations. 

V. INFRARED BEHAVIOR OF THE GLUON PROPAGATOR 

In Figures ||-|| we plot the data for the gluon propagator [see eqs. (||) and (^)] as a function of the lattice momentum 
p{k), defined in eq. (|l(]|), for different lattice volumes V and couplings /3. Our data confirm previous results |7|-[lC|] 
obtained in the strong-coupling regime for the four-dimensional case:M the gluon propagator is decreasing as p{k) 
decreases, provided that p{k) is smaller than a value pdec- Also, as in four dimensions, the lattice size at which this 
behavior for the gluon propagator starts to be observed increases with the coupling /?. In particular, in the strong- 
coupling regime, this propagator is clearly decreasing as p{k) goes to zero, even for relatively small lattice volumes 
(see Figures 2 and ph. On the contrary, for /3 > 3.4 (see Figure 0), this propagator is increasing (monotonically) in 
the infrared limit for V — 16'^, while it is decreasing at the largest lattice volume considered. 

Let us also notice that, at high momenta, there are very small finite-size effects, at all values of f3. The situation 
is completely different in the small-momenta sector, as already stressed above. Moreover, the value D{0) of the 
gluon propagator at zero momentum decreases monotonically as the lattice volume increases (see for example the 
case /3 = 5.0 in Figure^). These results suggest a finite value for D(0) in the infinite-volume limit, but it is not clear 
whether this value would be zero or a strictly positive constant. Therefore, the possibility of a zero value for D{0) in 
the infinite- volume limit is not ruled out. 

Finally, in Figure |g we plot (for j3 = 5.0 and 6.0 and V = 64^) the data for the gluon propagator D{k) with the 
choice k = (0,0, fc*), together with the data for k = (1, l,kt)- In both cases the two sets of data seem to fall on a 
single curve, i.e. we see no sign of breaking of rotational invariance. 

As said in Section ||, with our definition of the gluon field [see eqs. (g) and (^)] the quantity a~^A'' (x) approaches 

(1/2)5 [A^^°"*)]j'j(a;) in the continuum limit a ^ 0, where [A^^°"*)]^(x) is the (unrenormalizcd) conventional vector 

potential. In the same limit, aDik) approaches f;^_D'^'^°"*^(fc)/4, where £'('^°"*^(fc) is the unrenormalizcd continuum 
gluon propagator. [We recall that, in three dimensions, D^'^°"*^(fc) has mass dimension —2 and g^ has mass dimension 



In particular see Figure 1 in Reference 



1.] In Figure || we plot aD{k) (in GeV~ ) as a function of the momenta a~^p{k) (in GeV) for three different values 
of the coupling (3: 3.4,4.2 and 5.0. In all cases we consider the largest lattice volume available (see Table ||). The 
data show good scaling in the region where finite-size effects are negligible, i.e. in the limit of large momenta (see 
Figure |6|). |j On the contrary, the scaling is poorer in the infrared limit as expected. Nevertheless, we can see that 
the giuon propagator is decreasing for momenta p < Pdec, and that the value of pdec (in physical units) is practically 
/3-independent. From our data (see Figure ^p) we can set pdec ~ 350 MeV. Let us notice that pdec corresponds to the 
mass scale ^^/^ in a Gribov-like propagator. 

Finally, for the same set of data we consider (see Figure 0) the plot of the dimensionless product p{k) D{k) as a 
function of a~^p{k). Since p{k) D{k) is decreasing for a~^p{k) < 700 MeV, we can say that the gluon propagator is 
less singular than p~^ in the infrared limit, in agreement with Zwanziger's prediction. We notice that the turn-over 
value pto ~ 700 MeV is in good agreement with the result obtained recently in four dimensions for the SU{3) group 
(see Figure 8 in Reference Q). 

VI. CONCLUSIONS 

We think that our data for the gluon propagator are very interesting. The prediction p6|- [l9|j2l| ] of a propagator 
decreasing for momenta p < pdec is clearly verified numerically for several values of the coupling /?, ranging from the 
strong-coupling regime to the scaling region. Moreover, as in the four-dimensional case 0-^, it appears that the 
lattice size at which this behavior for the gluon propagator starts to be observed increases with the coupling. This 
requirement of large lattice volumes could explain why a decreasing gluon propagator has never been observed in the 
scaling region for the four-dimensional case [p]-[l0[. 

As mentioned above, our data in the strong-coupling regime and in three dimensions are in qualitative agreement 
with results obtained at small /3 in the four dimensional case Q-Q- This strongly suggests to us that a similar 
analogy will hold — in the limit of large lattice volumes — for couplings (3 in the scaling region, leading to an 
infrared-suppressed gluon propagator also in the four-dimensional case. 
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FIG. 1. Plot of the average plaquette (W^i.i) as a function of the couphng /3. For comparison we also plot the leading 
strong-coupling expansion /3/4 and weak-coupling expansion exp( — 1//3). Error bars are not visible. 
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FIG. 2. Plot of the lattice gluon propagator D{k) [see eqs. (pi) and (H)] as a function of the lattice momentum p{k) [see eq. 
([Lffl)] for lattice volumes V = 16"^ (□) and V = 32"^ (*), with k — (0, 0, kt), at /3 = 0. Error bars are one standard deviation. 
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FIG. 3. Plot of the lattice gluon propagator D{k) [see eqs. (@) and (H)] as a function of the lattice momentum p{k) [see 
eq. ^] for lattice volumes V = 16^ (□) and V = 32^ (*), with k = (0,0, fct), at: (a) /3 = 1.0; (b) /3 = 1.6; (c) /3 = 2.2; (d) 
/3 = 2.8. Error bars are one standard deviation. 
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FIG. 4. Plot of the lattice gluon propagator D{k) [see eqs. (B) and (0)] as a function of the lattice momentum p{k) [see eq. 
(0)] for lattice volumes V = 16^ (□), V = 16^x32 (+), V = 32^(*), V = 32^x64 (o), and V = 64^ (O), with k = (0, 0, fct), at: 
(a) /3 = 3.4; (b) /3 = 4.2; (c) /3 = 5.0; (d) /3 = 6.0. Error bars are one standard deviation. 
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FIG. 5. Plot of the lattice gluon propagator D{k) [see eqs. (fel) and (H)] as a function of the lattice momentum p{k) [see eq. 
@] for lattice volume V = 64^, with k = (0,0, kt) (*), and k = (1, 1, kt) (□), at: (a) /? = 5.0; (b) /3 = 6.0. Error bars are one 
standard deviation. 
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FIG. 6. Plot of the gluon propagator aD{k) in physical units (GeV ^) as a function of the momenta a ^ p{k) (in GeV) 
with k = (0,0, kt) and k = (1, l,kt), for /3 = 3.4 and V" = 32^ (□), /3 = 4.2 and V = 64^ (*), /3 = 5.0 and V = 64^ (O). In the 
second figure only the infrared region a~^ p{k) < IGeV is considered. Error bars are one standard deviation. 
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FIG. 7. Plot of the dimensionless product p{k) D{k) as a function of the momenta a ^ p{k) (in GeV) with k — (0, 0, kt) and 
k = (1, 1, kt), for f3 = 3.4 and V = 32^ (a), (3 = 4.2 and V = 64^ (*), f3 = 5.0 and V = 64^ (O). In the second figure only the 
infrared region a^^ p{k) < IGeV is considered. Error bars are one standard deviation. 
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TABLE I. The pairs {P,V) used for the simulations, the number of configurations, the number of HOR sweeps used for 
thermalization, the number of HOR sweeps between two consecutive configurations used for evaluating the gluon propagator, 
and the parameter pso used by the stochastic overrelaxation algorithm. 



/3 


V 


Configurations 


Thermalization 


Sweeps 


Pso 


0.0 


16-^ 


200 


4 


2 


0.85 


0.0 


32^ 


100 


4 


2 


0.90 


1.0 


le'' 


200 


1100 


100 


0.84 


1.0 


32^ 


100 


1100 


100 


0.93 


1.6 


16'' 


200 


1650 


150 


0.81 


1.6 


32^ 


100 


1650 


150 


0.83 


2.2 


16^ 


200 


2200 


200 


0.75 


2.2 


32''* 


100 


2200 


200 


0.76 


2.8 


16-^ 


200 


2750 


250 


0.72 


2.8 


32^ 


100 


2750 


250 


0.75 


3.4 


le'' 


200 


3025 


275 


0.69 


3.4 


32^ 


100 


3025 


275 


0.72 


4.2 


16'' 


200 


3300 


300 


0.66 


4.2 


32^ 


100 


3300 


300 


0.70 


4.2 


64^ 


56 


3300 


300 


0.72 



5.0 
5.0 
5.0 
5.0 
5.0 



16^ 
16^x32 

32^ 
32^x64 

64^ 



1000 
200 
170 
100 

54 



3575 
3575 
3325 
3575 
2275 



325 
325 
325 
325 

325 



0.63 
0.65 
0.68 
0.65 
0.69 



6.0 


16^ 


200 


3850 


350 


0.61 


6.0 


32^ 


150 


3350 


350 


0.67 


6.0 


64^' 


97 


3850 


350 


0.71 



TABLE II. For each coupling (3 we report the value of the average plaquette (Wi,i), together with the volume V and the 
number of HOR sweeps used for the analysis. Also, for (3 > 3.4, we report the tadpole-improved coupling /?/, the string tension 
^/a in lattice units, and the inverse lattice spacing a~^ in GeV. Error bars for the string tension and the inverse lattice spacing 
come from propagation of errors. Error bars for (VKi,i) are one standard deviation, evaluated taking into account the value of 
the integrated autocorrelation time Ti„t,Wi i for the plaquette, namely the variance is multiplied by 2Tint,Wi i (see eq. (3.7) in 
Reference pi). 







V 



Sweeps 



{Wi,i) 



Pi 



a-i (GeV) 



0.0 
1.0 
1.6 
2.2 
2.8 
3.4 
4.2 
5.0 
6.0 



32^ 
32^ 
32^ 
32^ 
32^ 
32^ 
64^ 
64^ 
64=* 



198 


0.000110(113) 








9900 


0.241650(16) 








14850 


0.373147(13) 








19000 


0.493302(12) 








24500 


0.595483(11) 








26950 


0.672720(9) 


2.28725(3) 


0.506(13) 


0.87(2) 


14149 


0.741862(4) 


3.11582(2) 


0.387(8) 


1.14(2) 


16164 


0.786877(3) 


3.93438(1) 


0.314(5) 


1.40(2) 


25135 


0.824783(3) 


4.94870(2) 


0.254(4) 


1.73(2) 



TABLE HI. For each coupling (3 > 3.4 we report the lattice spacing in fm, the largest lattice volume Vmax, the corresponding 
physical volume in fm^, and the smallest non-zero momentum that can be considered for that lattice (in GeV). 



/3 



a (fm) 



Vmax 



O, Vmax 



(fm^ 



^n (GeV) 



3.4 
4.2 
5.0 
6.0 



0.226(6) 
0.173(3) 
0.140(2) 
0.114(1) 



32-' 
64^ 
64^ 
64^ 



7.2'^ 
11. 1^ 
9.0^ 
7.30^ 



0.171 
0.112 
0.137 
0.170 
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